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Abstract: It is now well-known that vacuum polarisation in QED can lead to super- 
luminal low-frequency phase velocities for photons propagating in curved spacetimes. 
In a series of papers, we have shown that this quantum phenomenon is dispersive 
and have calculated the full frequency dependence of the refractive index, explaining 
in detail how causality is preserved and various familiar results from quantum field 
theory such as the Kramers-Kronig dispersion relation and the optical theorem are 
realised in curved spacetime. These results have been criticised in a recent paper by 
Akhoury and Dolgov [1], who assert that photon propagation is neither dispersive 
nor necessarily causal. In this note, we point out a series of errors in their work 
which have led to this false conclusion. 



1. One of the important simplifications in our analysis of photon propagation in 
curved spacetime is the insight that, in the eikonal and 'weak curvature' approxi- 
mations (see section 2), the background may be replaced by its Penrose plane wave 
limit around the null geodesic describing the classical trajectory. The first main 
claim of [1] is that such plane waves are too simple to manifest the vacuum polar- 
isation induced modifications to photon propagation discovered by Drummond and 
Hathrell [2]. This claim is simply not true: in fact plane waves have precisely the 
data that is encoded in the Drummond-Hathrell result. 

The origin of superluminal low-frequency propagation is the effective action for 
QED in curved spacetime [2]: 

£ = -^F^.F^'^ + [a.RF^^F^^- + a^R^^F^^F^ + a^^R^^^pF^^ F^^) +■■■ . 

(1) 

The dots indicate that this effective action is the first term in a derivative expan- 
sion, so results deduced from it are valid only for low-frequency propagation. The 
refractive index derived from (|I]) for a photon with wave-vector k'^ = uk^, where u 
is the frequency, is 

riij = 6ij + ^^{ciSijR^^M" + C2Rf,iujk^k'') , (2) 

where the indices i,j label the two spacelike polarisation directions. The constant 
coefficients ci,C2 are simply related to the known coefficients 02,03 in the effective 
action. If we now introduce a null coordinate u along the direction fc^, then 

nij = 6ij + —^^{ci^ijRuu + C2Ruiuj^ ■ (3) 

So the Drummond-Hathrell result depends on the curvature components Ruiuj and 

uu- 

Now, a plane wave metric in Brinkmann coordinates takes the form 

ds^ = 2du dv + hij {u)y' du^ + dy' dy' . (4) 

where {u, v) are null coordinates and are transverse spacelike coordinates. The non- 
vanishing components of the Riemann tensor describe the wave profile, Ruiuj = hij. 
So if we apply the Drummond-Hathrell formula (|^) to a wave propagating along 
the geodesic u with v = y^ = 0, that is k^ = u, k^ = and fc* = 0, we find 
a non-vanishing effect involving precisely the curvature components Ruiuj that are 
non-trivial in a plane wave spacetime. 

Of course, this is apparent in our work [3-6], since we have demonstrated that 
the low-frequency limit of our full, dispersive refractive index formula reproduces 
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the original Drummond-Hathrell result. The error in [1] is in their eq.(12), which is 
derived simply from the equation of motion in a Ricci flat plane wave background. 
This should read: 



with 0,3 = f ^aa, leading (since of course g™ = Quu = 0) to the usual result and 
not, as quoted in [1], 



The error is just a simple mistake of raising/ lowering indices with the off-diagonal 
metric in light cone coordinates u, v. However, it leads the authors of [1] to the 
important but manifestly wrong conclusion that the Penrose limit fails to capture 
the Drummond-Hathrell effect. 

The intuitive reason why the Penrose limit captures the essential geometry is 
most evident in the worldline, or proper-time, representations of the propagator 
where it is clear that the vacuum polarisation depends on the geometry of geodesic 
fluctuations around the original photon trajectory [4,5]. This is precisely what is 
captured by the Penrose limit and is encoded in the Van Vleck-Morette matrix, 
which plays a key role in our analysis. 

More recently [7], we have also made an independent check on the validity of 
the Penrose limit. In the case of a background spacetime dS^ x R, we can evaluate 
the refractive index, or spectral density, exactly. Then, taking the appropriate WKB 
limits (section 2), we recover precisely the result we obtain by using the Penrose limit 
directly, in this case a conformally-fiat symmetric plane wave with R^^ < 0. This 
gives a highly non-trivial confirmation that the Penrose limit does indeed retain the 
essential geometry of the full background field relevant for determining the vacuum 
polarisation corrections to photon propagation. 

2. Photon propagation with vacuum polarisation in curved spacetime is characterised 
by three length scales and it is important to be clear about the regimes of curvature 
and frequency where our results are valid. Let 1 / ^/TZ denote the generic curvature 
length scale (so TZ can represent a curvature or derivatives of curvature), u the photon 
frequency and m the electron mass. We work in the 'weak curvature' limit \/TZ <^ 
m (so the curvature scale is much greater than the electron Compton wavelength) 
and in the usual geometric optics, or 'eikonal' approximation u ^ \/TZ (so the 
photon wavelength is much less than the curvature scale). This defines the WKB 
limit in which the Penrose limit is a good approximation to the general background 
spacetime. 



2dud,Ak + {g^^ - 8d3Rukuk)dlAk = 0, 



(5) 



{2 + AdsR^kuk)dud,Ak + g^^dlAj, = 0, 



(6) 
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An important parameter in our analysis is the remaining dimensionless variable 
u\/TZ/w? and we need to consider all values of this to calculate the full frequency 
dependence of the refractive index. Causality is related to the limit n{oo) since this 
is identified with the wavefront velocity of the light wave. Since ref. [1] misrepresents 
the role of this parameter in our analysis of the high frequency behaviour of n{ijj), 
and therefore wrongly criticises the implications of our results for causality, we spell 
out the parameter dependence of our results here in explicit detail. 

Our result for the refractive index, taking scalar QED for simplicity, is (see, 
e.g. eq.(4.42) of [5]) 

a I , , / m 



with 

dt 







1- A{u,u-t)J det A (u, u-t) 



The curvature dependence is embedded in the Van Vleck-Morette matrix A(m, u — t) 
which is a function of \/TZt. For small t, it can be expanded in schematic form as 

OO 

A(«,«-t) ^ l + 5^7^"(M)^2" , (9) 

n=l 

where to linear order, the relevant curvature component is Ruiuj- (The vector nota- 
tion takes account of the polarisation dependence, i,j = 1,2, and will be dropped 
from now on for simplicity.) A simple rescaling t ^ t/z, or alternatively t — )■ t/y/R, 
shows immediately that eq.(|^ can be written in the (clearly equivalent) forms: 

„ = 1 + ^%f(^) = 1 + ^^^g(^) . (10) 

where the function F (or G) is calculated non-perturbatively. This allows us to access 
both the low and high frequency limits of the refractive index. 

At low frequencies, we can expand F as a power series in small ^^^^ using the 
VVM expansion (H). This starts with a term of giving the refractive index 

in the schematic form 

„ = 1 + -JL(i + of^)V (11) 



The first term, independent of frequency, reproduces the Drummond-Hathrell result. 
Of particular interest is the next term in the expansion which, if present, gives 
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a contribution to the imaginary part of the refractive index Im n{u) of the form 
uduR/m'^- These contributions can be reproduced by an extension of the DH effective 
action, discussed below. 

At high frequencies, we can use the second form of (0) and expand G for small 
^j^. Here, the first term is found to be of 0{1) and the expansion may contain 
logarithms as well as powers. We therefore find 

ayR. /-ir,\ 

n = 1 + + ... , (12) 

vr u 

with n(oo) = 1 as expected in a causal theory. 

All this is illustrated quite explicitly in a simple example discussed in [5], the 
conformally flat symmetric plane wave, where we have an exact analytic expression 
for ^{z) in @. In this model, the wave profile function is hij = cr'^dij, so \/TZ = a. 
We find: 

. .(£,(l,^)_£,e,^_l). (13) 



For large z i.e. low frequency, the digamma function has the expansion 

Z \ , Z a B2n /2cr\2". 

n=l 

where i?2n are the Bernoulli numbers, and so 



^(,) = _i^fi+on'). (15) 



3 2; V \z 
Substituting in (0) gives the refractive index [5] 

„H = i_^4(i + o(£!£!)). (16) 

907r V V / / 

Notice that in this case the expansion has no term linear in u and indeed Im n[u) = 
in this model. 

In the opposite limit of small z, i.e. high frequency, the relevant expansion of the 
digamma function is 

n=0 ^ ' 

giving 



m = a(-l--^log-+,j+0^- 
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The refractive index is therefore [5] 



n{u) 



1 - 



Unco 



a (7 




) 



) 



(19) 



where the first correction also includes the logarithm from (|l^). This example demon- 
strates precisely how both the low and high frequency limits of the refractive index 



refutes the claim of ref. [1] that the corrections to the low-frequency value n(0) van- 
ish in the high-frequency limit leaving n(oo) 7^ 1 with the associated problems with 
causality. 

3. In flat spacetime, an important role is played by the Kramers-Kronig dispersion 
relation 



In this case, Im n{u) > by virtue of the optical theorem, which relates it to the 
forward scattering cross-section, so n(oo) < ^(0) or equivalently fph(oo) > fph(O). 
A superluminal low-frequency phase velocity in flat spacetime would therefore imply 
the wavefront velocity fph(oo) exceeds c, violating causality.^ 

However, in curved spacetime neither the Kramers-Kronig dispersion relation 
nor the optical theorem remain true in their normal flat spacetime forms, and the 
usual insights based on (^) are simply wrong. Nevertheless, ref. [1] continues to use 
( pOD with the inevitable mistaken conclusions. Moreover, a confusion is made (eq.(2) 
of ref. [1]) between purely geometric and vacuum polarisation effects on Im n{u). 

To explain this second point, note that the eikonal approximation for the solution 
to the 0{a) corrected wave equation is (e.g. [6]) 



^An interesting, but quite separate, question is whether causahty would be violated in a model 
in which the strong equivalence principle is violated explicitly by tree-level curvature couplings and 
we can use a Drummond-Hathrell type action directly to analyse signal propagation and causality. 

In ref. [8], it was pointed out that this is related to 'stable causality' and depends on the global 
nature of the background spacetime. The question is whether the spacetime still admits a global 
Killing vector which is timelike with respect to the extended light cones of the DH action. Ref. [1] 
disputes this, repeating the proposal of Dolgov and Novikov [9] that the superluminal DH effect 
permits the construction of a 'time machine' which can be realised, for example, by two black holes 
in relative motion where each exhibits superluminal propagation in its vicinity. However, it has 
already been explained in [10], section 5, (see also [11]) that this proposed time machine does not 
work, precisely because of the fact that the DH effect is dependent on the local spacetime curvature 
(a fact dismissed in [1] as a 'trivial complication'), which invalidates the argument of [9]. 





(flat spacetime) . 



(20) 




(21) 
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where V is a. Rosen coordinate. The amphtude A{x) satisfies the classical equa- 
tion dulogA = —9 along the null geodesic, where 9 is the optical scalar in the 
Raychoudhuri equations [6] which describes the expansion or contraction of the null 
geodesic congruence describing the wave propagation. The amplitude can therefore 
increase or decrease due to this classical geometric effect, which should be clearly 
distinguished from the 0{a) change to the amplitude which would be induced by 
an imaginary part of n{u;u)) in (|2T|). This would correspond to genuine dispersion 
due to particle creation if Im n{uj) > or the more subtle photon 'undressing' effect 
described carefully in ref. [7] if Im n{uj) < 0. It is only this latter effect that is re- 
lated to vacuum polarisation, the Kramers- Kronig dispersion relation and the optical 
theorem. The assertion in [1] that such effects are 'clearly negligible' in this context 
completely misses the point, since the whole phenomenon of vacuum polarisation 
induced corrections to photon propagation takes place at this order. 

The form of the Kramers-Kronig dispersion relation which remains valid in 
curved spacetime is [5] 



This assumes only that n{u] u) is bounded at infinity and analytic in the upper half 
cj-plane, as required by causality. In order to recover the conventional form (|20|), 
we need to assume also that we have translation invariance along the geodesic, so 
that n{ijj) is an even function of u, and that it satisfies real analyticity, n{u*) = 
n{uj)*. This is described in full detail in [5]. One of the key discoveries of our 
work is that this second property does not hold in general, due to a novel analytic 
structure of the refractive index, or vacuum polarisation, in curved spacetime. This 
is due to additional singularities and branch cuts (evident in (pTSD) related to the 
existence of conjugate points on the null geodesies, i.e. points which may be joined 
by a continuous family of geodesies infinitesimally close to the original one. (In the 
world-line formalism, these correspond to zero-modes in the fiuctuations around the 
classical null geodesic path.) 

The novel analytic structure we have uncovered in curved spacetime may also 
give rise to a non-perturbative imaginary part for the refractive index, distinct from 
the uduR/m'^ type terms described above. These contributions to Im n{u) can arise 
even in spacetimes where the curvature is translation invariant along the photon's 
null geodesic, as e.g. the Ricci-fiat symmetric plane wave [5]. 




(22) 
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4. These entirely non-perturbative effects cannot be seen in a local effective action 
of the Drummond-Hathrell type, even when it is extended to include higher deriva- 
tives to all orders [6,12,13]. However, this is sufficient to capture the perturbative 
imaginary parts proportional to uduR/m^- Since this is also relevant for the central 
argument of ref. [1] (see section 5), we recall some details of the derivation here. The 
generalised effective action is: 

+ + t,R^,F^^F\ + ^^R^^y^pF^^F^p) 

+ -^^{a[RD^F^^^D,F\ + h[R,,D^F^^^D,F'^'' 



+ b^2RuuD^F'^D''F^, + t,R,,D^^D^F^.FP^ + ^^R.^^.D^F'^pD^F^^ 



(23) 



In this formula, the o^, 6„, are known 'form factor' functions of three operators, 
i.e. 



r)2 r)2 r)2 

-> _ ^ (^W ^(2) ^(3) 






(24) 



where the first entry D'^^^ acts on the first following term (the curvature), etc. 

The refractive index derived from (p3D is [13] 

, , ^ ^ a 1 {2iuk-D\^ a 1 /2iuk ■ D\ ^ 
nij{uj) = dij + dij -Ci{ — ]Ruu H ^ — ]Ruiuj , y^^^) 

where the constant coefficients Ci, C2 of (^ are replaced by functions of the operator 
k ■ D du, which describes the variation of the curvature tensors along the classical 
null geodesic. The linear terms in the expansion of ci and C2 give the contributions 
to Im n{u) of C(^-;^;^) discussed following (0). We have checked that these agree 
with those derived using our VVM-based method [6]. 

To understand the origin of this expression, consider as an example the term 

J d^x^R,,,,F'^^D'F'P (26) 



in the effective action, incorporating a single power of D'^ from the form factor. The 
corresponding contribution to the equation of motion for D^F^'^ is 



1 



-2{D^R^\p)D^F^P - A{D„R^\p)D''D^F^P - AR^\,D^ D'^ F^p 



(27) 
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The key point is that only the second term survives at leading order in the eikonal 
approximation, where it contributes to the C2 term in (|25|) . The first term is sup- 
pressed by a power of IZ/w? while the third is of 0(/c^), which is suppressed by a 
power of a since the photon is on-shell before the radiative corrections. By contrast, 
the second term is of only of relative order Oi^^^^^), the now familiar parameter. 

5. The main constructive part of ref. [1] is based on a 1994 paper by Khriplovich 
[14] which argues on the basis of a proposed general form for the graviton-photon- 
photon vertex that all corrections to the refractive index beyond the low-frequency 
contribution n(0) vanish. Although at first sight plausible, it is apparent that this 
argument misses the contributions of the form k^^D^R/rn^ ~ uduR/m^ that we 
have found in both the effective action and VVM approaches. This was pointed out 
already in ref. [13]. Of course, Khriplovich's method could not in any case detect the 
non-perturbative contributions following from the occurrence of conjugate points on 
the photon's null geodesic. 

Although it is obvious from our explicit construction of n{u) in many examples 
that Khriplovich's argument cannot be correct, the loophole in the analysis of [14] is 
in fact quite subtle. Following [14], we define the graviton-photon-photon vertex in 
momentum space on a flat background: 

{27cy5^^\q + k + k')T{q,k,k') = W''{q) f dz e''-''{t{k)\T^,{z)\t{k')) , (28) 



where e{k) is the photon polarisation and h^v{q) is the external graviton field. Then, 
according to [14] (see also [2,15]), the most general form of the vertex to 0{l/m^) 
can be written as 

T{qXk') = a,{q\k',k")R\q)F,,{k)F^^{k') + a^iq', k\ k")R';^^{q)F^\k)F\{k') 
+ a,{q\k\k'')R%^^{q)F^^{k)F^^{k') , 

(29) 

where F^u{k) = k^ey{k) — k^e^ik) and R^{q) is the Fourier transform of the Ricci 
scalar evaluated with the metric identified with the graviton field, etc. The aj(g^, /c^, fc'^) 
are form factors. The similarity with the effective action (EBI) is clear. Higher-order 



Lorentz structures of 0{l/m^) with four momenta can be read off from (p3D (see 
also, e.g. [16]). 

The argument of [1, 14] now proceeds by claiming that the imaginary part of 
r, with the photons taken on-shell, only has a contribution from = so that 
Imr(g^, 0, 0) ~ ^{q^)- Since the vertex has no dependence on the photon frequencies, 
it follows that the refractive index itself must be independent of frequency and that 
n(uj) = n{0) for all u, i.e. the Drummond-Hathrell effect is non- dispersive. 
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To see what is wrong with applying this argument in curved spacetime, notice 
that for a slowly varying background gravitational field, the Fourier transform R{q) 
of the curvature is in fact very singular at g = 0. If we expand (dropping indices for 
clarity) : 

R{z) = R{(}) + zf'df.RiO) + ... , (30) 

then 

R{q) = R(0)S^^\q) - ^^,R{0)^6^'\q) + ■■■ . (31) 

To see the effect of this singular behaviour involving derivatives of the delta function, 
substitute (0) and integrate a typical vertex term over q: 

J d^q a{q\ k\ {k + qf) R{q)F{k)F{-k - q) 

= a{0,k\k^) R{0)F{k)F{-k) - a'^^\0,k^,k^) 2ik''di,R{0) F{k)F{-k) (32) 
+ a{0,k^k^) td^R{0) F{k)^F{-k) , 

where a^^-* denotes differentiation w.r.t. the third argument. Of these terms, the 
first is obviously just the zero-momentum Drummond-Hathrell contribution, while 
the third is lower order in the eikonal approximation. The second term, however, 
can contribute at leading order even on-shell and gives rise to an imaginary part for 
the refractive index proportional to k^df^R{0), i.e. precisely the frequency- dependent 
contribution of uduR/m^ type that we have already identified. It arises in essentially 
the same way as illustrated for the effective action in (|26| ) , (p^ ) . 

What this shows is that when the Khriplovich argument is applied to curved 
spacetime backgrounds, it is necessary to keep careful track of the full singularity 
structure, including derivatives of the momentum-space delta function associated 
with the graviton insertion. Yet again, we find that a careful analysis shows that 
the refractive index for photon propagation on curved spacetime is dispersive, with 
a low-frequency expansion (0) which may have an imaginary part proportional to 
k^'D^R/m^. 

6. Finally, the question of whether QED in curved spacetime is causal is ulti- 
mately determined by whether the commutator, or Pauli- Jordan, Green function 
iG^i,{x,x') = {0\[A^{x), A^{x')]\0) vanishes outside the light cone. This was demon- 
strated in ref. [5] by an explicit construction, contradicting the claims in [1]. In fact, 
we were able to show that the commutator in scalar QED in the Penrose plane wave 
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limit, incorporating vacuum polarisation, can be written as 




(33) 



X Aij {u, u-t) \/detA({i,{i - t) G 




; x I + • • • , 



whee the omitted terms are independent of the curvature. Here, G {m^\x,x') is the 
commutator Green function for a free massive scalar field, which clearly has support 
only on or inside the light cone. Since Gij{x,x') therefore vanishes outside the light 
cone, causality is manifest. 

It is interesting to relate this to the refractive index. If we take x — {u,V,0,0) 
and x' — {u' — > — oo, 0, 0, 0) to be two points on the classical photon trajectory, 
we find (note that since we are taking u > u', the commutator and retarded Green 
functions coincide): 



The properties we have established for the refractive index, in particular analyticity 
in the upper half a;-plane, now show that the commutator (retarded) Green function 
vanishes when V <0, i.e. outside the light cone. 

In summary, despite the many erroneous claims to the contrary in ref. [1], we 
have clearly established [3-6] that photon propagation in QED in curved spacetime 
is indeed dispersive and causal. 
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